tic;
% Objective
% f=[-3,-5,-3];
% f=zeros(3,1);
f=randi([2001],1000,1);

% Equality Constraints
% A=[1,-1,1;2,3,6];
% b=[15;30];
% A=zeros(2,3);
A=randi([2001],1000);
b=randi([2001],1000,1);
% Equality Constraints
% Aeq=[2,1,0];
% beq=[30];
% Aeq=zeros(1,3);
Aeq=randi([2001],1,1000);

% beq=zeros(1,1);
beq=round(rand(1)*2001);

% lower bound
lb=zeros(1000,1);

% upper bound
% ub=[15;10;5]; 
% ub=Inf(3,1);      
% ub(1)=15;
% ub(2)=10;
% ub(3)=5;
ub=[];
toc;

% solver
tic;
[x,fval,exitflag,output]=linprog(f,A,b,Aeq,beq)
% [x,fval,exitflag,output, lambda]=linprog(f,A,b,Aeq,beq,lb,ub)
toc;
